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The phase behaviour of amphiphilic multiblock copolymers with a large number of blocks in semidilute solutions is studied 
by lattice Monte Carlo simulations. The influence on the resulting structures of the concentration, the solvent quality and the 
ratio of hydrophobic to hydrophilic monomers in the chains has been assessed explicitely. Several distinct regimes are put in 
evidence. For poorly substituted (mainly hydrophilic) copolymers formation of micelles is observed, either isolated or connected 
by the hydrophilic moieties, depending on concentration and chain length. For more highly substituted chains larger tubular 
hydrophobic structures appear which, at higher concentration, join to form extended hydrophobic cores. For both substitution 
ratios gelation is observed, but with a very different gel network structure. For the poorly substituted chains the gel consists of 
micelles cross-linked by hydrophilic blocks whereas for the highly substituted copolymers the extended hydrophobic cores form 
the gelling network. The interplay between gelation and phase separation clearly appears in the phase diagram. In particular, 
for poorly substituted copolymers and in a narrow concentration range, we observe a sol-gel transition followed by an inverse 
gel-sol transition when increasing the interaction energy. The simulation results are discussed in the context of the experimentally 
observed phase properties of methylcellulose, a hydrophobically substituted polysaccharide. 
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1 Introduction 

Amphiphilic block copolymers 1 are particularly versatile 
macromolecules because they allow for a rich variety of dif- 
ferent structures. They may be built from many different kinds 
of monomers with different chemical properties. Their length 
and the number of blocks of each species can be tuned at 
will, from di- and triblock to multiblock copolymers -. Their 
architectures can be linear, branched or star-like, the blocks 
may be distributed randomly or regularly. Also, copolymers 
may interact through hydrogen bonding, hydrophobic interac- 
tions or electrostatic forces, thus resulting in different, pos- 
sibly competing, types of behaviour in solution, depending 
on the physico-chemical parameters such as temperature, pH, 
ionic strength and type of solvent 2 " 4 . Due to their dual affinity 
for polar and apolar solvents these macromolecules can self- 
assemble in order to minimize their free energy. This gives 
rise to a wide range of structures and properties, from swollen 
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copolymers in a good solvent to microphase formation, gela- 
tion and phase separation in a poor solvent. Amphiphilic block 
copolymers have attracted a great deal of attention from the 
biomedical field 5 , with applications such as drug delivery vec- 
tors^, nanoparticle stabilizers, nanoreservoirs, emulsion stabi- 
lizers, wetting agents, rheology modifiers 7 8 or as injectable 
scaffold materials for tissue engineering 9 . 

In the present work our focus is on the specific case of 
strictly alternating linear multiblock copolymers consisting of 
two kinds of monomers, hydrophobic (H) and polar (P). This 
may be viewed as a generic, simplified model for the under- 
standing of the structural aspects and the phase behaviour of 
hydrophobically substituted polysaccharides 10 , an important 
example of which is methylcellulose. Methylcellulose is ob- 
tained by methylation of some hydroxyl groups of the cellu- 
lose polymer. It results in alternating blocks of more or les s 
hydrophobic monomers, whose proportions and lengths de- 
pend on the degree of methylation of each glucose cycle (from 
0, hydrophilic to 3, hydrophobic). There is experimental evi- 
dence that the substituted monomers are indeed arranged in a 
block pattern—"—. Hence we expect that the simulation of a 
hydrophobic/polar multiblock copolymer can give insight into 
the structure and properties of methylcellulose as a function of 
thermodynamic condition s 14 i 15 . 

There is a vast body of theoretical and experimental studies 
on diblock and triblock copolymers- 7 *^"—. Much fewer stud- 
ies deal with the simulation of amphiphilic multiblock copoly- 
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mers, and most of them address the questions of their struc- 
ture and kinetics in dilute solutions where only intramolecu- 
lar interactions play a role 22,23 . It was convincingly demon- 
strated that multiblock copolymers not only form intramolec- 
ular chains of micelles 24,25 but also tubular structures and 
intramolecular bilayers— . The formation of this variety of 
intramolecular structures is mainly influenced by two pa- 
rameters, the effective interaction strength between the H 
monomers and the hydrophobic substitution ratio. It is a con- 
sequence of the balance between the entropic penalty of con- 
straining the polar P monomers and the energetic gain due to 
the self-assembly of the hydrophobic H monomers. The for- 
mation of intramolecular pearl-necklaces of micelles in multi- 
block copolymers was justified theoretically by Halperin 27 . 
For the case of amphiphilic graft copolymers, the existence of 
intramolecular pearl-necklaces of micelles and cylindrical mi- 
celles was argued theoretically and directly observed in sim- 
ulation studies 28,29 . Also, it was shown experimentally that 
amphiphilic graft copolymers in dilute solutions form either 
unimolecular rodlike chains of micelles or unimolecular rod- 
like micelles, depending on the quality of the solvent 30 ' 31 . 

Regarding semidilute solutions, where both intra- and inter- 
molecular interactions are present, Anderson and Travesset 32 
performed simulations of the phase behaviour of pentablock 
copolymers and found different regimes as a function of con- 
centration and solvent quality of the hydrophilic blocks: a 
swollen gel, giant micelles, lamellar and lamellar-perforated 
phases were observed. Glotzer et al— showed in a two- 
dimensional model calculation that multiblock copolymers, 
unlike diblock copolymers, micro-phase separate and form a 
gel network. Gindy et al. 34 focussed on the role of the number 
m and length n of the blocks in (H„P„)„, multiblock copoly- 
mers containing the same number of H and P monomers. They 
showed that these copolymers may either evolve by micro- 
phase separation or by macroscopic phase separation, depend- 
ing on the m/n ratio. From the experimental side microphase 
separation of graft copolymers was studied as a function of 
the hydrophilic and hydrophobic block lengths 35 . From a The- 
ological point of view Regalado et al. 36 advanced evidence 
for three viscoelastic regimes in semidilute solutions of multi- 
block copolymers, depending on concentration, hydrophobic 
ratio and hydrophobic block length. A study of the impact of 
pH and temperature on the supramolecular self-assembly of 
pentablock copolymers in solution 37,38 showed that increas- 
ing the pH leads to an increase in the hydrophobicity which, 
in turn, triggers the formation of spherical micelles which 
eventually become cylindrical (also called wormlike). Lin- 
tuvuori and Wilson 39 focussed on the behaviour of rod-coil 
multiblock copolymer melts when cooling a system through 
the order disorder transition (ODT). They gave evidence for 
the formation of micellar, nematic, lamellar and gyroid phases 
and showed that the length of the rods influences the tempera- 



ture of the ODT. Recently, the influence of a block profile that 
varies along the copolymer has been studied. The resulting 
micelles significantly deviate from the ones obtained for di- 
block copolymers 40 . The dynamics of block copolymer sys- 
tems has been critically discussed by Denkova et al.—, paying 
particular attention to the consequences of slow dynamics on 
experiments and simulations. 

In the present work we focus on the behaviour of linear 
multiblock copolymers in semidilute solutions, where both 
intra- and intermolecular interactions come into play. In con- 
trast with previous investigations, in our model each copoly- 
mer is made of a large number of hydrophobic and hydrophilic 
blocks (up to 15 of each). 

The role of the polymer concentration and the strength of 
the hydrophobic effect on the phase behaviour is studied. Fur- 
thermore, the role of the substitution ratio is investigated. Dif- 
ferent patterns are identified (including spherical and cylin- 
drical micelles, gel formation and phase separation) and the 
transitions between these different structures are assessed. 

The paper is organized as follows: Section |2]describes the 
simulation method, the model and the simulation details. The 
properties computed for gaining an understanding of the be- 
haviour of the system are defined in Sec. [3] The results are 
presented in Sec. [4] for two representative substitution ratios. 
They are discussed and compared with previous experimental 
and simulation studies in Sec. [5] 



2 Simulation model and method 
2.1 Simulation method 

In the present work an efficient lattice MC model is used, 
which we introduced in detail in a previous article 26 . Polymer 
chains are discretized on a face-centered cubic lattice, thus al- 
lowing for a large number of nearest neighbour sites and bond 
angles. Periodic boundary conditions are imposed. A lattice 
site may be occupied by zero, one or two monomers. Ad- 
jacent monomers in a chain occupy the same site or nearest 
neighbour sites on the lattice (double occupancy is restricted 
to neighbouring monomers along a polymer, see Hugouvieux 
et al— for full details). Thus a lattice site represents a 
volume of two monomers. The excluded- volume constraint 
is strictly respected by all monomers and the polymers are 
never allowed to cross each other. Dynamics is implemented 
through attempted nearest neighbour hops of randomly se- 
lected monomers. Moves are accepted if the final configura- 
tion again respects bond length and excluded volume restric- 
tions. For the interacting hydrophobic monomers a trial move 
furthermore has to satisfy the Metropolis criterion to be ac- 
cepted 26 . 
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2.2 Model 

We consider semidilute solutions of linear, regularly alternat- 
ing multiblock copolymers made of two kinds of monomers, 
hydrophobic (H) or hydrophilic/polar (P). The solvent is not 
dealt with explicitly, and the hydrophobic effect is mimicked 
by an effective nearest-neighbour attractive interaction be- 
tween the H monomers, which implicitely expells the polar 
solvent from the vicinity of the hydrophobic monomers. This 
attractive force acts between any two H monomers that oc- 
cupy the same or two nearest-neighbour sites (called interact- 
ing monomers in the following) and is given by the interaction 
energy E, < 0, expressed in units of the thermal energy kgT . 
Increasing the absolute value of corresponds to a strength- 
ening of the H-H interaction and hence is equivalent to a low- 
ering of the temperature (except for the case of systems with 
a lower critical solution temperature (LCST), for which it cor- 
responds to an increase in temperature 42 ). No difference is 
made between intra- and intermolecular interactions between 
H monomers. The P monomers neither interact explicitely 
with the H or the P monomers nor with the solvent. The pat- 
tern of the blocks in the chains is described as {Hb h Pb p )„, 
with Bh the length of a hydrophobic block, Bp the length of 
a polar block and n the number of times such a diblock is 
repeated in the chain. In this study n takes values between 
5 (10 blocks) and 15 (30 blocks). Hence the focus is on sys- 
tems of copolymers with a large number of blocks. The degree 
of polymerization (total number of monomers) of the chain is 
N,„ =nx (Bn + Bp). The hydrophobic substitution ratio is de- 
fined as P su y = B H" Bp ■ Distances (wave vectors) are expressed 
in units of lattice spacing (2k over lattice spacing). 

2.3 Simulation details 

Simulations were performed for different values of concentra- 
tion, solvent quality, substitution ratio and copolymer chain 
length. The concentration (j) is defined as the fraction of 
the total available volume which is occupied by monomers. 
Because of the double occupancy the total volume V = L? 
is twice the number of available lattice sites N s = jL 3 , 
where L is the periodicity. Simulations were performed 
on monodisperse many-chain systems with chain lengths 
N„, = 100,300, for different concentrations (j) ranging from 
0.005 to 0.24, and for different interaction energies £, from 
-0.3 to -0.6. This enables a broad range of structural features 
to be encompassed, from dilute good solvent conditions to 
concentrated solutions of strongly interacting copolymers. 
Also, different substitution ratios P su b, chain lengths N m 
and numbers of polymer N p give insight into the phase 
behaviour and finite-size effects. The simulations are per- 
formed in the canonical NVT ensemble, with fixed number 
of monomers N, volume V, and temperature T, For each 
set of conditions, several independent simulation runs are 



performed (typically 8 to 16), starting from distinct random 
initial configurations, generated by a random polymerization 
process. Time is measured in Monte Carlo step (MCS) where 
one MCS corresponds to N m x N p trial moves. The systems 
are equilibrated for some 10 7 MCS. The relaxation of the 
time auto-correlation function of the end-to-end vector of the 
polymers is used as a criterion for equilibration. During the 
production run configurations are analysed periodically (every 
2. 10 s to 10 6 MCS depending on Ej and 0) in order to compute 
the properties of the system. Typical production runs range 
from 10 7 to 8.10 7 MCS, depending on the relaxation time 
of the system. Simulations were performed for two values 
of the substitution ratio, P su b = 0.2 and 0.5, representative 
of different copolymer architectures. The simulation details 
are summarized in Table Q] the concentration being varied 
through changes of the system size. 

Table 1 Parameters used in the simulations 



Psub B H Bp N p N m (j) 

0.2 4 16 60 300 0.005 -0.04 

0.2 4 16 100 100 0.005 -0.04 

"05 5 5 50 fiX) 0.005 - 0.05 

0.5 5 5 300 100 0.005 -0.24 



3 Properties 

Direct visualization of the configurations gives valuable qual- 
itive information about the structure of the system. The calcu- 
lation of experimentally measurable quantities from the simu- 
lations allows for a quantitative assessment of these properties 
and of the influence of the different physical parameters. The 
calculated quantities are defined in the following. 

3.1 Clusters 

A cluster is defined as a set of interacting copolymers. Two 
copolymers interact whenever at least one H monomer of the 
first copolymer is nearest neighbour to an H monomer of the 
second copolymer. A cluster comprises all the monomers (H 
and P) of the involved chains. If polymer Pi interacts with 
P2 and, P3 interacts with P2, all three polymers belong to 
the same cluster. Thus connectivity through the hydrophilic 
blocks within a copolymer is part of the cluster definition. 

3.2 Hydrophobic cores 

A hydrophobic core is defined as a set of interacting H 
monomers, ignoring the P monomers and the connectivity of 
the chains. If the H monomer Hi interacts with H2 and H2 
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interacts with H3 all three belong to the same core. Note that a 
cluster may contain several hydrophobic cores. If we consider 
a single chain, it always belongs to a single cluster for connec- 
tivity reasons, but it may contain as many hydrophobic cores 
as the number of H blocks in a chain, especially in good sol- 
vent conditions. We use the following notation to define dif- 
ferent useful observables: A^ is the number of H monomers in 
a given core, A^, is the number of chains whose H monomers 
are part of a given core and N^' 1 is the number of H monomers 
of a given polymer that belong to a given core. All the distri- 
butions W(N) presented below are normalized as follows: if 
N occurs 0(N) times, W(N) = N.O(N)/Zn> N'.0(N') (e.g. 
W(Nfj) is the fraction of H monomers in cores of size N^,). 



qualitative differences of the resulting structures for different 
substitution rates. Therefore we will present data for two dis- 
tinct substitution rates, each being representative of one sce- 
nario. We distinguish between a low substitution rate case, 
where the screening of the hydrophilic coronae is complete 
and a high substitution rate case where the hydrophilic screen- 
ing is incomplete. In general the properties of a copolymer 
do not only depend on the substitution ratio, but they depend 
explicitely on both the hydrophilic and hydrophobic block 
lengths in a non-trivial way. We have considered different 
possibilities studying a single copolymer 2 '', here we choose 
one case for each scenario and call it the low and the high 
substitution ratio case, respectively. 



3.3 Percolation 



4.1 Low substitution ratio, P SU b = 0.2 



A cluster is considered to be percolating if it spans the whole 
simulation box, i.e. if one or more of its dimensions are 
larger than the size of the box. A configuration percolates 
if it contains at least one percolating cluster. For a given set 
of conditions the probability of percolation, F 'percolation . is de- 
fined as the fraction of percolating configurations. We take 
Ppercoiution = 0.5 to be the onset of gelation. 

3.4 Structure factor 

Calculating the scattering intensity of the simulated sys- 
tems yields quantitative information on structural correla- 
tions which may be compared directly with experimental data 
from scattering experiments—. Because the structure of the 
monomers is neglected in our simulations, there is equivalence 
between the scattering intensity and the structure factor and 
the two terms will be used interchangeably. The structure fac- 
tor is defined as: 



N N 

s{q) = E E ( ex pH<7- (n-h)}) 

;=U=l 



(1) 



where N is the number of monomers in the system, q a vec- 
tor of the reciprocal space and q its norm, f> and the posi- 
tions of monomers j and k, respectively. S(q) gives informa- 
tion about the structural correlations in the system on different 
length scales. It is useful to compute partial structure factors 
defined as the separate contributions to S(q) of only the H or 
only the P monomers. This is eventually to be compared with 
experiments using labelled chains-^. 

4 Results 

It has been shown for isolated copolymers, that the sub- 
stitution rate determines the pattern developping, through a 
balance between energetic (hydrophobic) and entropic (hy- 
drophilic) contributions 26 . As will be shown below, there are 



We first present data for a low substitution ratio with an excess 
of hydrophilic monomers. The influence of the concentration 
and the interaction energy on the formation of hydrophobic 
cores were monitored separately. Figures [T] and |2] show the 
hydrophobic core size distribution either as a function of con- 
centration or interaction energy. The core sizes are expressed 
in number of H monomers per core, A^, which is necessarily 
a multiple of the number of H monomers in a block, Bh- The 
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Fig. 1 Size distribution of the hydrophobic cores, W(Nfj), for 
different concentrations <|> at E t = —0.51, P su b = 0.2, B# = 4, 
Bp = 16, N p = 100, N m = 100. 

core size distribution as a function of concentration is shown 
in Figure Q] and has the typical characteristics of a system of 
surfactant micelle s 44 ' 45 . In the present context a micelle is a 
structure with an H core surrounded by a P shell. Its shape 
is either spherical, tubular or lamellar. The appearence of a 
bump in the size distribution of the hydrophobic cores Nfj 
(when the slope at the inflexion point becomes positive) in- 
dicates the distinct presence of a well-defined core size, the 
micelle core. This is used as a definition for the critical mi- 
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celle concentration (CMC) 46 . A similar criterion has been 
used for monitoring the formation of surfactant micelles—. In 
FigureQ], at E t = -0.51, the CMC is (j> = 0.01 and the dis- 
tribution of core sizes becomes approximately concentration- 
independent above (]) = 0.02, with a stable peak in the mi- 
celle size distribution developping around 70 H monomers. 
Since each copolymer contains 20 H monomers, typical mi- 
celles necessarily contain hydrophobic blocks from different 
polymers, they are multimolecular. A somewhat different be- 
haviour is observed when increasing the interaction strength 
(see Fig. [2]) at a fixed concentration (j) = 0.02. There is a sys- 
tematic increase of the number of large micelles with increas- 
ing interaction strength while the maximum of the distribu- 
tion lies around 65 H monomers and increases slightly with 
energy. The number of isolated H blocks (i.e. H blocks that 
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Fig. 2 Size distribution of the hydrophobic cores, W{Nji), for 
different interaction energies at § = 0.02, P su i, = 0.2, Bh = 4, 
Bp = 16, N p = 100, N m = 100. 

do not interact with any other H block) decreases both with 
concentration and interaction strength. At fixed = —0.51, 
it decreases by a factor of w 1 .7 upon increasing the concen- 
tration from (j) = 0.005 to 0.035, and at fixed = 0.02, it even 
decreases by a factor of w 2.8 when increasing the interaction 
strength from £, = —0.47 to —0.55. Similar simulations have 
been performed for longer chains, N m = 300 and N p = 60 (see 
Table[T]i. The size distribution of the hydrophobic cores for the 
longer copolymers is essentially the same both with regard to 
its general shape and the value of the most probable core size. 
Consistent with the observation that the miceilar structure is 
solely stabilised by the balance between the energy of its hy- 
drophobic core and the entropy of the surrounding hydrophilic 
shell and does not depend on chain length, the multimolecular 
aspect of the core structure for N m — 100 at low concentration 
becomes essentially monomolecular for N m = 300. For the 
given set of parameters N m = 300 happens to be just the limit 
of the chain length beyond which one expects that a pearl- 



necklace of micelles starts to form 26 . 
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Fig. 3 Size distribution of the hydrophobic cores, W(Njj), for 
different concentrations at a stronger interaction energy Ej = —0.55, 
P sub = 0.2, B H = 4,B P = 16, N p = 100, N m = 100. Micellization 
occurs at any concentration, the CMC is equal to zero. 

In Figure [3] we show the evolution of the core size distribu- 
tion at P su b = 0.2 for a stronger interaction energy Ej = —0.55, 
for different concentrations. By contrast with £, = —0.51 mi- 
celles are now present at any concentration and the core size 
distribution becomes almost concentration-independent, with 
a peak between 60 and 70 H monomers. Slight deviations are 
visible at the lowest concentrations. Hence the difference be- 
tween the core size distributions at Ej = —0.5 1 and £/ = — 0.55 
is that in the latter case the CMC is equal to zero. This is 
because, at £, = —0.55, the effective concentration of the H 
monomers along a single copolymer is high enough such that 
they always collapse into a miceilar structure. Simulations 
with N p = 60 and N m = 300 give very similar results, with the 
notable exception that the miceilar peak is much more narrow 
at low concentration. With increasing concentration the peak 
becomes significantly broader (not shown). We attribute this 
to the fact that the commensurability between the number of 
H monomers in a copolymer and in a miceilar core is better at 
low concentration. In fact, when increasing concentration, the 
number of different chains participating in a core increases and 
the distribution develops a broad maximum at a value close to 
3 copolymers per core. Correspondingly, the distribution of 
N^' L , the number of H monomers in a core belonging to a given 
polymer, changes its character completely, as shown in Figure 
|4] At low concentration it is sharply peaked at Ng = 60, 
which is just the value when all H monomers of a copolymer 
belong to the same core. At high concentration there is instead 
a broad maximum around A^' c = 20. As the size of a core 
does not change when concentration is increased, the number 
of different cores connected by the hydrophilic P blocks of 
a polymer must increase with concentration. This is directly 
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Fig. 4 Size distribution of the subset of H monomers of a given 
copolymer belonging to a given core, W(N^' C ), for different 
concentrations at Ej = -0.6, P sub = 0.2, B H = 4,B P = 16, N p = 60, 
N m = 300. The maximum value on abscissa corresponds to the total 
number of H monomers in a copolymer. 



confirmed by the measurements of the radius of gyration of 
a copolymer, which increases from R\ = 19 at = 0.005 to 
flg = 44 at = 0.041. 

The phase diagram in Figure [5] summarizes the simulations 
performed at P su b = 0.2, as a function of concentration and 
interaction energy. Depending on the parameter values the 
system may form spherical micelles and/or a gel. The (£,-,0) 
conditions where the system contains micelles are determined 
using the inflexion point criterion of the core size distribution 
(blue line in Fig. |5). Geometric percolation is taken as the 
hallmark of gelation 47 (red line inFig.|5}. Four regimes can be 
distinguished: no micelles and no percolation; only micelles; 
only percolation; both micelles and percolation. Let us first 
consider the evolution of the system with respect to micelle 
formation. For Ej < —0.52 micelles are observed no matter 
how low the concentration. On the other hand, for Ej > —0.52, 
the formation of micelles depends on the concentration of the 
system as shown in Fig.Q] the lower the interaction strength, 
the larger the concentration needed for the onset of micelle 
formation (CMC) 46 . This relationship is almost linear. 

A snapshot of a configuration corresponding to the point A 
in Fig. |5]is shown at the top of Fig. [6] It shows a rather dilute, 
strongly fluctuating set of hydrophobic cores surrounded by P 
monomers, the micelles. Percolation sets in for concentrations 
larger than = 0.025. The shape of the sol-gel phase bound- 
ary is curved such that, in a narrow concentration range, the 
system first undergoes a sol-gel and then a gel-sol transitions 
when increasing the interaction strength at constant concen- 
tration. This can be understood in view of the results pre- 
sented in Fig. [3] Increasing the interaction strength increases 



the core size. Because of the compactness of the hydrophobic 
cores, this leads to an effective shortening of the chain which, 
in turn, prevents gelation. There is a large domain in which 
micelle formation and percolation overlap. Point B in Fig. [5] 
lies in this domain of Ej and values. Typical configurations 
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Fig. 5 Phase diagram of all systems simulated with short 
hydrophobic and long hydrophilic blocks: formation of micelles 
and/or percolation as a function of energy E, and concentration for 
P sub = 0.2, B H =4,B P = 16, N p = 100, N m = 100 and box sizes 
ranging from L — 62 to L = 126. Dots denote the parameter values 
where simulations were performed. Lines separate the different 
phases. Formation of micelles is observed in the blue area while 
gelation occurs in the red area. Note that around ~ 0.026 there is 
successively a sol-gel and then a gel-sol transition upon increasing 
\Ej\. In Fig. [6] snapshots of typical configurations are displayed, they 
correspond to points A and B in Fig. [3], 

in this zone show a set of micelles consisting of a hydrophobic 
core surrounded by and connected to each other with blocks of 
hydrophilic monomers (see Fig. [6] bottom). Complementary 
results (not shown) obtained for N p = 60 and N m — 300 have 
the same phase structure but with the percolation line shifted 
to lower concentrations. This is consistent with the fact that 
increasing the chain length facilitates cross-linking. 

Structure factors yield valuable information on the spatial 
organisation of the system ^ For copolymers it is particu- 
larly useful to calculate partial structure factors, as they give 
insight into the spatial correlations of each of the two kinds of 
monomers separately. This is shown in Fig. |7]for P su f, = 0.2 
for two different concentrations 0. The partial structure fac- 
tors are defined by equation Q] where the summations are 
now restricted to the P or the H monomers, respectively. At 
= 0.02 the strongly interacting system (E{ = —0.55) shows 
the typical features of spherical core-shell structures 49 ' 5Q with 
a Guinier plateau for the H monomers corresponding to the 
size of the cores and a lack of intensity for the P monomers at 
the same length scale. Also, a Porod representation of 5//(q r ) 
(not shown) exhibits a plateau corresponding to the presence 
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monomers shows more statistical noise at low q and the de- 
crease in intensity at intermediate q for the P monomers is 
less pronounced. The Porod representation of Sn(q) no longer 
shows the plateau characteristic of a sharp segregation of H 
and P monomers. Concomitant gelation and core formation in 
this system can explain these observations: the cross-linking 
of the micelles impedes the formation of well-defined core- 
shell structures with low interfacial tension. The increase of 
the large scale statistical fluctuations may be due to the more 
slowly relaxing gelling systems which implies that fewer inde- 
pendent configurations are being sampled. The partial struc- 
ture factors computed for longer copolymers, N p = 60 and 
N m = 300, confirm that at equilibrium the core size and shape 
do not depend on chain length. 
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Fig. 7 Partial structure factors S(q) for P s „i, = 0.2. Top: (|> = 0.02, 
bottom: (]) = 0.032. Left: P monomers only, right: H monomers 
only. At Ej = —0.55, the presence of cores is clearly visible in the H 
structure factors. 



Fig. 6 Snapshots of typical configurations for P su i, = 0.2, i.e. short 
hydrophobic and long hydrophilic blocks. Top: micelle phase, point 
A in Fig.[5](after 10 7 MCS). Bottom: gelled micelle phase, point B 
in Fig.|5](after 10 7 MCS). Note that both simulations are performed 
with the same numbers of monomers but at different concentrations, 
which explains the difference in box size. Red beads: H monomers. 
Black dots: P monomers. Some monomers seem isolated because of 
the periodic boundary conditions. 

of a sharp interface between the H and P monomers. This 
contrasts the same analysis for the more weakly interacting 
system at = —0.47 where no micellar core is visible. 

At higher concentration, (j) = 0.032, the structure is less 
well-defined: the Guinier plateau at £, = —0.55 for the H 



4.2 High substitution ratio, P su b = 0.5 

Probing systems with a higher substitution ratio allows a dif- 
ferent micellar regime to be explored. While for Bh = 4 and 
Bp = 16 well defined spherical micelles are stable in the en- 
ergy range considered, for Bh = Bp — 5 they are not. Figures 
[8] and [9] show the time evolution of the size distribution of 
the hydrophobic cores for Bh = Bp = 5. Data is presented 
for strongly interacting H monomers, deep inside the binodal 
domain (£,■ = —0.5), and for two concentrations, one below 
and one above the percolation threshold ((j) = 0.01 and 0.085). 
The size distributions at (]) = 0.01 and 0.085 look very differ- 
ent, nevertheless they have some common features, notably 
the presence of well defined peaks. These peaks appear for 
monomer numbers that are multiples of fifty, which is the to- 
tal number of H monomers in a polymer. This indicates that 
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most chains collapse to have all their H monomers in a sin- 
gle core, forming either a monomolecular micelle (size = 50 
monomers) or multimolecular micelles (multiples of 50). 
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Fig. 9 Same as Fig.[8]for = 0.085. A series of peaks 
corresponding to exact multiples of the number of H monomers in a 
polymer develops with time and is most pronounced for 
intermediate core sizes. 



Fig. 8 Time evolution of the size distribution of the hydrophobic 
cores, W[Njj), for high substitution rate, P su i, = 0.5, for <|> = 0.01, 
Ei = -0.5, P sub = 0.5, B H =B P = 5, N p = 300, N m = 100. The 
different curves correspond to different simulation times. Note the 
pronounced peak structure corresponding to exact multiples of the 
total number of H monomers in a copolymer. 

In the low concentration case (see Fig. |8j most cores have, 
at all times, a size which peaks sharply for multiples of fifty. 
Furthermore, the core sizes markedly evolve towards larger 
values with time. This is visible in the figure and is also 
confirmed from the distribution of N H ' , the number of H 
monomers in a given core and a given polymer (not shown): 
at all times the monomers of a given polymer belong to the 
same core and the cores evolve by aggregation of collapsed 
polymers. The initial collapse of a copolymer into a globule 
takes place on a much shorter time scale than the coarsening 
of the cores. 

At higher concentration, (j) = 0.085, at early times the core 
sizes are not discretized in units of fifty, indicating that the 
initial phase separation is no longer due to the collapse of in- 
dividual copolymers (Fig. But with increasing simulation 
time not only larger cores form as shown by the overall shift 
of the distribution towards higher values, but also peaks start 
to appear progressively. The peaks first appear between 10 7 
MCS and 2. 10 7 MCS and then shift to larger multiples of fifty. 
This is again related to the appearance of clusters of fully col- 
lapsed polymers : starting from 2.10 7 MCS, the peaks of the 
distribution for intermediate core sizes, corresponding to a few 
isolated copolymers, develop. This is the same behaviour as 
for <j) = 0.01. By contrast, for larger core sizes the peaks tend 
to be more smeared out. These results suggest that the initial 
core formation is constrained by the connectivity of the chains, 
then there is partial reorganisation to form copolymer globules 



and finally these globules further evolve into larger structures. 
For a very dense system with (j) = 0.24 at E, = —0.5, there is 
still a broad initial maximum and an evolution towards larger 
micelles, but at no time any sharp peaks appear (not shown). 

The way how the distribution of core sizes is affected by 
concentration is summarized in Fig.[l0] With increasing con- 
centration we observe a systematic decrease of the intensity of 
the peaks and a shift towards larger core sizes. For = 0.01 
the largest core contains up to 500 H monomers (10 chains 
or more) while for = 0.24 the largest cores may contain up 
to 10000 H monomers (200 chains or more). The distribu- 
tion of the number of chains participating in a core, W(N^ h ), 
shows a significant shift towards a larger number of chains as 
well as a broadening when increasing the concentration from 
= 0.01 to 0.24. For = 0.24 the character of the distribu- 
tion differs qualitatively from the lower concentrations: there 
appears a second maximum at very large core sizes typical of 
the presence of a percolating core structure. Visual inspection 
of typical configurations confirm the existence of a spanning 
network of branched cylindrical micelle cores. 

Fig. ITTI shows that the time evolution of the distribution at 
low and at high concentration is totally different. For = 0.01 
a well-defined distribution develops with a maximum around 
3 chains per core. For = 0.24 a very broad peak around 
eight copolymers per core develops early on and systemati- 
cally spreads out to larger values at later times while the num- 
ber of single H blocks and small copolymer cores does not 
change much. The most significant feature of these curves 
is the appearance beyond 2.10 7 of a secondary peak at large 
sizes, confirming the appearence of a percolating core struc- 
ture, as already found in Fig. [10] 

The coarsening of the hydrophobic cores is also visible in 
the evolution of the average core size (not shown). Initially, 
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Fig. 10 Size distribution of the hydrophobic cores, W(N£j), for 
different concentrations at E\ = —0.5, P m h = 0.5, Bjj = Bp = 5, 
N p = 300, N m = 100. 
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Fig. 11 Time evolution of the distribution of the number of different 
copolymers in a core, W(N£ h ), for low ((]) = 0.01, dashed lines) and 
for high ((]) = 0.24, full lines) concentration. Both distributions still 
evolve. Note the appearance of a second peak for large sizes at 
d> = 0.24. 



there is a rapid saturation at the core size corresponding to 
isolated spherical micelles. As the relaxation is not complete 
a gradual change corresponding to the reorganisation of the 
micelles is visible at all concentrations. The onset of this slow 
coarsening increases with increasing concentration. A similar 
relaxation behaviour is observed when monitoring the evolu- 
tion of the average of Ng, the number of H monomers in a 
polymer that belong to a given core. 

The phase diagram for the substitution rate P m b = 0.5 is pre- 
sented in Fig.[T2]for concentrations ranging from (]) = 0.005 to 
0.24 and energies ranging from —0.3 to —0.5. Let us first 
consider the upper part of the diagram where no micelles 
are formed. It comprises a range of concentrations above 



(j) = 0.02 to 0.03 where the system percolates. Percolation 
being defined geometrically and the interaction strength be- 
ing weak, the gelation criterion may not be meaningful in 
this domain. Below the line of micelle formation the percola- 
tion threshold separates strongly interacting non-gelling sys- 
tems at low concentrations from gels at high concentrations. 
The gelation threshold increases with increasing interaction 
strength. Snapshots of typical conformations of the system, 
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Fig. 12 Phase diagram for P sub = 0.5, B H = 5,B P = 5. The two 
data sets are N p = 50, N m = 100 with L ranging from L = 46 to 
L = 100 (circles) and N p = 300, N m = 100 with L between L = 50 
and L = 180 (squares) (same notation as in Fig. [5}- The points for 
different box sizes coincide, showing that finite size effects are 
negligible. Snapshots of typical configurations corresponding to 
points C and D are shown in Fig. Q~3] 



at the strongly interacting points C and D in the diagram, are 
shown in Fig. [13] At low (j), point C, core-shell micelle-like 
structures are formed, typically spherical or elongated and oc- 
casionally connected by hydrophilic blocks. At high concen- 
tration, point D, extended linear and/or branched tubular core 
morphologies are observed. They can be assimilated to worm- 
like micelles^ 1 . Two system sizes with a different number of 
chains (A^, = 300 and N p = 50) were compared to check for 
finite size effects. The phase diagram determined from both 
systems is the same which confirms that the results are not 
affected by system size. 

The signature of the structural features are also visible in 
the partial structure factors presented in Fig. [14] At low con- 
centration ((j) = 0.01) Su(q) behaves, for all q, in a similar 
way as Sh(<1) at P su b — 0.2 and (j) = 0.01. This is because in 
both cases the H monomers form spherical micelle cores. At 
(|) = 0.085 the shape of Sn(q) is also similar to the P sl ,i, = 0.2 
case for large q. For low q values it exhibits a q~ 1 57 behaviour 
instead of the Guinier shape observed for the (]) = 0.01 cases 
for P m b = 0.2 and P m i, = 0.5. This suggests that on a large 
scale the H cores effectively scale like a polymer in good sol- 



1-03 1 9 




**** ' 



a result of the ongoing association and reorganization of stable 
collapsed copolymers. 
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Fig. 14 Partial structure factors S(q) for P sub = 0.5. Top: (|> = 0.01, 
bottom: (]) = 0.085. Left: P monomers only, right: H monomers 
only. 




Fig. 13 Typical conformations for P m i, = 0.5. Top: simulation of 
point C in Fig.[[2] Bottom: same for point D, both after 6.10 7 MCS 
(same notation as in Fig.|6j. 



vent (this is not strictly true as the H cores are more rigid and 
ocasionally branch). The hydrophilic partial structure factor 
Sp(q) is, both for low and for high concentration, in agree- 
ment with the expected shell structure, with the minor differ- 
ence that for (j) = 0.085 the profile is flatter for small q values. 
The time evolution of the partial structure factors (not shown) 
indicates that on intermediate q scales the final conformation 
(the core-shell structure) develops rapidly. Aggregation then 
proceeds more slowly on large length scales (low q values), as 



5 Discussion 

The purpose of this study is the investigation of the different 
scenarios of semidilute solutions of amphiphilic linear multi- 
block copolymers. For simplicity regularly alternating hy- 
drophobic and hydrophilic monomers were considered. De- 
pending on the hydrophobic substitution ratio P su b, the con- 
centration (]) and the strength E, of the attractive force between 
the hydrophobic monomers, different types of structures and 
phases are found, including good solvent behaviour, micelle 
formation, gelation and phase separation, as summarized in 
the phase diagrams (see Figs. [5]and[T2]for different substitu- 
tion ratios). 

For low concentration and weak interaction between the H 
monomers (white area in the diagrams) the system forms nei- 
ther micelles nor a gel for any substitution ratio. With increas- 
ing concentration there is a sol-gel transition that is no differ- 
ent from gelation of homopolymers. 

When increasing the interaction strength while staying at 
low concentrations self-assembly takes place (the area be- 
low the blue line in the phase diagrams). A common feature 
of the self-assembled structures is the formation of cores of 
hydrophobic monomers surrounded by shells of hydrophilic 
monomers, as a result of the short-range attractive interac- 
tion between the H monomers. Core formation of the copoly- 
mers may be triggered either by increasing the interaction en- 
ergy or the concentration. The threshold for structure for- 
mation depends on the hydrophobic substitution rate and on 
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the block lengths, as expected from dilute solutions —: the 
lower the entropic contribution of the P monomers (the shorter 
the P blocks in the chain), the lower the attractive interac- 
tion strength needed for triggering the association of the H 
monomers. 

Depending on the hydrophobic substitution rate and the 
chain length, copolymers show different types of core- 
shell structures, ranging from spherical micelles and pearl- 
necklaces of spherical micelles to linear and branched tubu- 
lar micelles. It is influenced by the P layers surrounding the 
H cores and the connectivity between the H and P moieties. 
The protection of a hydrophobic core depends directly on the 
lengths of the H and P blocks: Decreasing the number of hy- 
drophilic monomers in the chain leads to a thinner P shell as 
shown by the more pronounced intermediate peak in S(q) of 
the P monomers (compare Figures |7]and [14}. 

For thick P shells (P sl ,b = 0.2 in our study), the phase separa- 
tion of the H monomers into cores is stopped by the protective 
P layer as can be seen in Fig. [3] An increase of the core size 
would require P monomers to be included in the core, which 
is energetically unfavorable, or the entropic penalty from den- 
sifying the shell of P monomers would be too strong. Con- 
versely, for large values of P su b (P S ub = 0.5 in our study), the 
rather thin P shell does not protect fully against further asso- 
ciation of the H blocks. This lack of protection leads to elon- 
gated core-shell structures, elliptical at low concentrations and 
tubular and/or branched at higher concentration. Put differ- 
ently, the thickness of the outer P layer determines which mi- 
cellar structure stabilizes. A result similar to our simulations 
has been found experimentally 37,38 where spherical core-shell 
structures become elongated and eventually precipitate when 
increasing the hydrophobic! ty of the copolymers. 

At P si ,b = 0.2, when the micelles are well protected, the 
phase structure is quite similar to micellization in surfactant or 
diblock copolymer solutions, where, in a narrow range of 
increasing concentration induces micelle formation (CMC) 46 . 
In all simulations the micelles have a well-defined size dis- 
tribution which barely evolves with increasing concentration, 
while the number of micelles increases linearly. Also, the 
higher the concentration (in the range of Ej between -0.47 and 
-0.52), the lower the interaction strength needed for micelle 
formation. This is not surprising in the light of the fact that, for 
the chain lengths and simulation parameters considered, stable 
micelles are made of several copolymers. There are, however, 
noticeable differences with respect to surfactants. For strong 
interaction energies, copolymer micelles form at any concen- 
tration, the CMC is zero. This is because intramolecular mi- 
celles form at any polymer dilution as the connectivity locally 
forces a finite concentration of H monomers. For the same 
reason, the micelle phase boundary (see Fig.|5]and Fig.fl4l is 
rather flat. 

A parallel can be made between the effects of £, and P su b- 



Increasing the interaction between H monomers at low values 
of P su b where pearl-necklaces of micelles are present leads to 
the formation of tubular structures with the same radius. In- 
deed, for a given radius imposed by the block structure, the 
tube is energetically more favorable because of its lower sur- 
face to volume ratio. Upon increasing the interaction strength, 
this energetic gain eventually overcomes the entropic penalty 
of accommodating the P monomers around a tube. On the 
other hand, starting from a tubular structure and increasing 
the length of the P blocks (decreasing the ratio of hydropho- 
bic monomers in the chain) at constant Ej would lead to the 
break-up of the tube into several cross-linked spherical mi- 
celles, a pearl-necklace. 

The connectivity of the P and H moieties has a strong in- 
fluence on the core-shell structures. For low substitution rates 
(Psub = 0.2) the length of the H blocks determines the radius 
of the H cores, in the case of highly substituted copolymers 
(Psub = 0.5), the length of the H blocks limits the diameter 
of the tubular micelles. This relationship has been shown to 
hold for long copolymers in dilute solutions 26 . It also works 
quite well for more concentrated solutions: for P su b = 0.2 and 
Ei = —0.55 the core diameter estimated from S(q) for the H 
monomers is 3.8, to be compared with the block length 4; for 
Psub = 0.5 and = —0.5 the tube diameter is 5.2, for a block 
length of 5. 

The present simulations explore concentration effects due 
to aggregation of collapsed globular copolymers. At low sub- 
stitution rates (P sl ,b = 0.2) and low concentrations, the self- 
assembled structures are shown to be either multimolecu- 
lar ((H+PuOs, N m = 100) or approximately monomolecular 
((//4Pi6)i5, N m = 300). From the comparison of the two cases 
one concludes that the connectivity has little influence on the 
micelle stabilisation. This result is consistent with the findings 
of Gindy et al, whose simulations showed that increasing the 
number of diblock motifs in copolymer chains leads to a de- 
crease in the most probable number of chains involved in a 
core- 3 - 4 -. It was also shown experimentally by Zhou et al— 
that for a constant number of diblock motifs in the chain, 
increasing the length of both blocks while keeping the hy- 
drophilic/hydrophobic ratio constant induces a transition from 
multimolecular micelles to monomolecular micelles. There- 
fore the number of H monomers needed for the formation 
of stable core-shell structures is determined primarily by the 
energy-entropy balance and not by the chain length or the 
block structure. 

At low concentration and for the thermodynamic parame- 
ters and the chain lengths considered in our simulations, the H 
blocks of a copolymer all collapse into a single micelle core. 
The structural changes upon increasing concentration is very 
different for low and for high P su b- For P sl ,b = 0.2 the core- 
shell structure does not change significantly when increasing (]) 
but the H blocks of a given copolymer are now shared between 
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several cores thus providing hydrophilic links between differ- 
ent cores (this is particularly prominent for N m = 300 where 
the cores are monomolecular at low (j) and intermolecular at 
high (j)). At P su b = 0.5 the core structure is not preserved when 
increasing (j) because the low concentration spherical/elliptical 
micellar cores aggregate as a whole to form larger cores. 

For gelation we use the static percolation threshold as the 
criterion. For weak energies, i.e. \Ej\ < 0.48 for P su j, = 0.2 
and \Ej\ < 0.35 for P su b = 0.5, when no micelles are formed, 
percolation occurs for purely geometric reasons and does not 
necessarily lead to a physical gelation transition. For energies 
where phase separation of the H monomers sets in, percolation 
does correspond to gelation. At P su b = 0.2 gelation is observed 
either by increasing concentration or upon strengthening of the 
interactions. A particularity of the concentration driven transi- 
tion is that gelation not only consists of connecting nearby mi- 
celles, but also because increasing the concentration induces 
the hydrophobic blocks of a copolymer to spread over several 
micelle cores which effectively increases size and entangle- 
ment of the copolymers. The gel looks rather weak as it con- 
sists of hydrophilic cross-links between spherical hydrophobic 
cores. Interestingly, in a narrow concentration range the sys- 
tem undergoes a reverse gel-sol transition at quite strong in- 
teractions. This transition is due to the copolymer contraction 
upon strengthening the interaction, as directly visible from the 
chain radii which decreases from 7? 2 w 23 at Ej = —0.54 to 

ps 17 at Ej = —0.6. For P sl ,t, = 0.5 the percolating network 
is very different: The spanning cluster is made of a contin- 
uous structure of hydrophobic monomers, covered by a shell 
of hydrophilic monomers. The evolution of the system from 
spherical micelles to tubular structures to a spanning network 
proceeds by coalescence since at high density the P shell can- 
not prevent the H monomers of different cores from interact- 
ing. Furthermore, the tubes tend to be straight because the P 
shell at the tube ends is less dense which favors growth at the 
ends. Also, at P sll t, = 0.5 no reverse transition from gel to sol 
occurs at strong energies because increasing the energy does 
not significantly affect the topology of the hydrophobic tubu- 
lar cores. Consistently with the different cases for P si ,b = 0.2 
and 0.5, simulations of pentablock copolymers in solution by 
Anderson and Travesset 32 have provided evidence for what 
they call a swollen gel (a network of loosely connected mi- 
celle cores) and for a giant cylindrical micelle (a single tubu- 
lar core made of several copolymers) when increasing the ratio 
of hydrophobic to hydrophilic monomers in the chain. Also, 
Gindy et al.~ observed in simulations, at P su y = 0.5, the for- 
mation of a different kind of gel network for what they call 
microscopically or macroscopically phase separating copoly- 
mer, depending on the lengths of the H blocks, the P blocks 
and the copolymers. 

The present simulations distinguish, depending on the pa- 
rameter values, between two classes of micelles, well pro- 



tected ones {P S ub = 0.2, spherical) and poorly protected ones 
(Psub = 0.5, tubular). For the chain lengths considered each 
micelle consists of one or a few fully collapsed copolymers. 
Our low concentration calculations 26 for much longer chains 
(N m > 300) show that under well protected conditions pearl 
necklace structures and under poorly protected conditions ex- 
tended tubular (wormlike) structures appear. Extensive cal- 
culations for these cases at high concentration are computa- 
tionally demanding but one may guess the gel structure that 
develops in this case upon increasing concentration. For pearl 
necklaces occasional bonding of micellar pearls, by the same 
mechanism as the hydrophilic bonding of spherical micelles in 
our P su b = 0.2 simulations in Fig.|6l will lead to a loosely and 
hydrophilic ally connected network of necklace chains. For the 
tubular case the aggregation process is no different from the 
Psub = 0.5 one in Fig.[13]but with a coarser network structure 
because with increasing chain length the elementary micelles 
simply are more elongated. 

Regarding the phase behaviour of methylcellulose, the 
present model gives a possible explanation for the two differ- 
ent types of gels observed experimentally: a weak, transpar- 
ent gel at low temparature and a strong, turbid gel at higher 
temperature. The clear gel at relatively low concentrations 
and weak interaction strength then corresponds to the mi- 
cellar phase connected by hydrophilic strands and the tur- 
bid gel at higher temperature (because of the LCST, higher 
temperature corresponds to stronger interactions) corresponds 
to a gel formed by a network of thick hydrophobic strands. 
In an experimental situation, and particularly for methylcel- 
lulose, the distinction between hydrophobic (strongly inter- 
acting) and hydrophilic (non-interacting) blocks is of course 
not as clearcut as in the model, therefore the hydrophilic 
monomers will also contribute to the structure formation. For 
methylcellulose the degree of methylation of each glucose 
cycle ranges from (hydrophilic) to 3, (hydrophobic)—"—. 
Therefore, in the strong gel regime intermediate degrees of 
substitution provide some interaction between the strongly in- 
teracting hydrophobic cores. The strand formation observed in 
cryo-TEM images suggests that rather stiff bundles form in the 
strongly phase separating part of the phase diagram 53 ' 54 . Thus 
the multiblock model can explain the main features of methyl- 
cellulose phase diagram. The observed bundling 53 could be 
included by introducing a weak attraction between the hy- 
drophilic monomers in our model. There is one puzzling ex- 
perimental feature that is not reproduced in the simulations: 
in the weak gel region there is a gel to sol transition upon in- 
creasing concentration 14 . Possible explanations for this coun- 
terintuitive observation could be a more complicated interac- 
tion structure (e.g. involving hydrogen bonds) or that the The- 
ologically determined gel-sol transition is not a macroscopic 
phenomenon but the effect of intermediate scale structural re- 
organisations. 
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6 Conclusions 

In the present study a generic lattice copolymer model is 
used to assess the phase behaviour of amphiphilic multiblock 
copolymers in semidilute solution, as a function of the solvent 
quality of the H monomers, the length of the hydrophobic and 
hydrophilic blocks in the chain, and the copolymer concentra- 
tion. The simulations provide evidence for two distinct types 
of gel structures, depending on the stability of the micelles 
that form when the hydrophobic monomers phase separate. 
For chains of spherical micelles, also called pearl necklaces, a 
weak, hydrophilic ally linked gel forms when, with increasing 
concentration, the micelle cores exchange hydrophobic blocks 
belonging to different copolymers. Such a gel is a direct con- 
sequence of the multiblock nature of the copolymer. For un- 
stable micelles a strong gel builds up when spherical or tubular 
elementary micelle cores join to form larger linear or branc hed 
structures. By lowering the solvent quality one passes from 
the weak gel to the strong gel phase. This provides a possible 
explanation for the weak and the strong gel phase observed in 
methylcellulose. 
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